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ABSTRACT 


The  Effective  Drag  Functions  derived  by  the  Naval  Weapons  Laboratory  for  the 
computation  of  gun  ballistic  data  in  digital  gunfire  control  systems  are  presented. 
The  accuracy  and  derivation  of  the  functions  are  given  as  well  as  the  least-squares 
coefficients  computed  for  3-  and  5-inch  projectiles.  FORTRAN  listings  of  the  three 
programs  used  in  the  computation  of  the  coefficients  are  given  in  the  appendices. 
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I.  INTRODUCTION 


Reference  1  states:  “A  gun  weapon  system  is  functionally  composed  of  the 
following  elements:  (1)  a  fire  control  system,  (2)  one  or  more  gun  mounts 
organized  as  a  battery,  and  (3)  ammunition  supply  arrangements.  With  guns,  fire 
control  can  be  considered  as  the  practical  application  of  exterior  ballistics  to  ensure 
that  the  projectile  will  hit  the  target.  In  any  weapon  system,  the  fire  control 
problem  is  essentially  that  of  getting  the  weapon  or  projectile  to  hit  the  target  and 
explode,  or  at  least  to  explode  when  the  target  and  the  weapon  or  projectile 
approach  close  enough  for  maximum  damaging  effect.  In  gun  weapon  systems,  the 
problem  is  to  lay  (aim)  the  gun  in  such  a  way  that  the  projectile  will  hit  or 
approach  close  to  the  target.  Unlike  other  types  of  weapon  systems,  once  the 
projectile  has  left  the  gun  muzzle,  nothing  further  can  be  done  to  affect  its  course. 
The  gun’s  fire  control  problem  ideally  must  be  solved  before  the  gun  fires.  To  solve 
the  gun  fire  control  problem,  it  is  necessary  to  deal  with  three  main  types  of 
variables:  (1)  exterior  ballistics,  (2)  target  position  and  relative  motion  of  target  and 
own  ship,  and  (3)  inherent  corrections  necessitated  by  the  physical  characteristics  of 
the  weapon  system  and  the  medium  on  which  its  carrier  floats.” 

As  indicated  above,  one  of  the  major  roles  of  the  fire  control  system  is  to 
accurately  predict  the  exterior  ballistic  performance  of  the  projectile  since  nothing 
can  be  done  to  affect  its  course  once  it  has  exited  the  muzzle.  Ultimately  then,  the 
accuracy  with  which  the  system  can  predict  the  projectile’s  trajectory  provides  an 
upper  limit  on  the  capability  of  the  system  to  achieve  its  designed  purpose, 
i.e.,  incur  maximum  damaging  effect  on  the  target.  From  the  instant  the  projectile 
leaves  the  muzzle  of  the  gun  until  it  ends  its  flight  by  impact  or  explosion,  the 
projectile  trajectory  is  affected  by  the  following  factors:  (1)  momentum,  (2)  gravity, 
(3)  air  resistance,  (4)  wind,  (5)  drift,  (6)  earth’s  rotation  (Coriolis),  and  (7)  earth’s 
curvature.  Thus,  the  fire  control  system  must  make  provision  for  these  factors  in 
aiming  the  gun  such  that  the  projectile  will  intersect  the  path  of  the  target. 

Normally,  a  range  table1  (i.e.,  firing  table)  is  constructed  for  each  projectile 
based  on  exterior  ballistics  theory  and  information  gained  from  the  analysis  of 
instrumented  flight  tests  of  the  projectile.  The  range  table  is  based  on  a  set  of 
nominal  or  standard  firing  conditions  for  the  weapon  and  provides  differential  effects 
for  variables  such  as  initial  velocity,  projectile  weight,  air  density,  air  temperature, 
wind,  arid  gun  and  target  motion.  The  table  gives  elevation  angle,  time-of-flight,  and 
drift  for  standard  firing  conditions;  using  the  differential  effects,  corrections  can  be 
applied  to  the  data  for  standard  firing  conditions  to  obtain  a  gun  elevation  angle 


^ec  Reference  2,  3,  and  4. 


and  an  azimuth  angle  for  any  combination  of  firing  conditions  that  may  exist.  It 

should  be  noted  that  the  range  table  is  an  approximation  to  the  trajectory  of  the 

projectile.  Solutions  for  non-standard  firing  conditions  have  inherent  errors  in  them 
due  to  the  non-linearity  of  the  differential  effects  given  in  the  table  and  the 
interaction  of  multiple  effects  operating  simultaneously.2 

Reference  1  defines  fire  control  computers  or  rangekeepers  as  electronic  or 
electromechanical  devices  into  which  are  fed  the  mathematical  variables  in  the  fire 
control  problem,  and  which  yield  solutions  in  the  form  of  control  settings  (e.g. ,  gun 

orders)  required  for  the  weapon  and  launcher  to  have  maximum  effect  on  the 

target. 

All  operational  gunfire  control  computers  in  the  U.  S.  Navy  are  analog  and  use 
either  mechanical  cams  and/or  potentiometers  and  other  electronic  circuitry  to 
essentially  reproduce  the  minimally  essential  columns  of  the  gun  range  table.3  The 
later  versions  are  electrical  analog  using  modern  solid  state  circuitry  to  approximate 
relatively  simple  polynomial  functions  which  regenerate  portions  of  the  range  table. 
The  functions  are  necessarily  simple  and  cither  neglect  some  of  the  variables  which 
affect  projectile  flight  or  only  correct  for  their  variations  in  an  approximate  fashion. 
These  simplifications  are  required  on  the  basis  of  economy  and  space  limitations. 
For  example,  none  of  the  analog  systems  provide  inherent  corrections  for  Coriolis  or 
consider  the  effects  due  to  the  interactions  of  several  var  ations  from  standard  firing 
conditions  occurring  at  once.  It  should  be  noted  that  under  certain  conditions 
provision  is  made  to  correct  for  factors  which  have  smaller  influence  on  the  flight 
of  the  projectile  than  those  that  are  neglected.  Interaction  effects  under  certain 
conditions  can  influence  the  flight  of  the  5-inch  MARK  41  projectile  as  much  as 
5  mils,  a  value  generally  considered  larger  than  the  dispersion  of  the  round  itself. 

By  the  early  1960’s,  computer  technology  had  advanced  to  the  state  that 
digital  computers  were  being  built  that  could  be  used  to  solve  the  tactical  fire 
control  problem.  The  substitution  of  digital  for  analog  computers  showed  promise  of 
much  greater  flexibility  and  sophistication  and  therefore  greater  accuracy  in  the 
solution  of  the  gunfire  control  problem.  For  example,  in  the  area  of  exterior 
ballistics  it  permits  a  much  more  accurate  prediction  of  the  trajectory  of  the 
projectile  than  was  possible  with  the  analog  computer.  Ideally,  the  best  approach  to 
predicting  the  trajectory  of  the  weapon  in  the  fire  control  system  is  to  numerically 
integrate  the  differential  equations  of  motion  within  the  fire  control  computer. 
However,  this  is  an  iterative  time-consuming  process  which  has  not  yet  been  proven 
practical  in  the  gunfire  control  computer,  especially  in  light  of  the  high  frequency 


o 

Sec  Reference  5. 
3Sec  Reference  6. 


of  solutions  needed  to  achieve  system  stability.  The  most  feasible  alternative  is  to 
develop  ballistic  functions  to  approximate  the  trajectory  data  obtained  from 
numerical  integration  of  the  equations  of  motion.  With  the  digital  computer,  this 
can  be  done  with  a  degree  of  accuracy  which  is  practically  impossible  with  the 

analog  computer. 

Since  its  inception,  the  Naval  Weapons  Laboratory  (NWL)  has  been  the  Navy’s 
prime  facility  for  exterior  ballistics  for  surface-launched  unguided  weapons,  including 
guns.  In  the  early  1960’s,  NWL  was  asked  by  the  Naval  Ordnance  Systems 
Command  (NOSC)  to  serve  as  technical  advisor  to  Lockheed  Electronics  Company 
(LEC)  in  the  development  of  the  U.  S.  Navy’s  first  digital  gunfire  control  system. 
Range  tables  and  specialized  ballistic  data  were  provided  to  Lockheed  in  support  of 
their  efforts  which  included  implementing  the  range  tables  in  the  system.  The 

ballistic  functions  used  in  the  MARK  86  GFCS  are  partitioned  cubic  least-squares 
polynomials  fit  to  the  essential  columns  of  both  the  surface  and  anti-aircraft  range 
tables  for  5"/54  caliber  gun.  The  ballistic  solution  is  essentially  one  of  applying 
differential  corrections  to  range  and  time  of  flight  to  obtain  gun  orders  such  that 
the  projectile  will  strike  the  target  under  any  reasonable  non-standard  firing 

conditions  that  may  exist.  The  solution  is  analogous  to  that  described  in  the 
introductory  information  to  the  surface  range  table. 

Initially,  the  MARK  86  GFCS  was  designed  as  a  surface  only  system;  therefore, 
no  provision  was  made  to  include  ballistic  solutions  for  anti-aircraft  (AA)  targets, 
that  is,  targets  at  non-zero  position  angles.  When  the  decision  was  made  to  add  an 
AA  capability,  surface-type  range  tables  were  fit  with  polynomials  (identical  to  the 
procedure  for  the  surface  table)  for  various  position  angles.  To  obtain  an  AA 

solution  with  the  system,  surface-type  solutions  are  obtained  for  three  position 
angles  which  bracket  the  target  position  angle;  quadratic  interpolation  is  then 
performed  among  the  three  surface-type  solutions  for  the  position  angle  of  the 
target. 

This  procedure  requires  a  large  amount  of  computer  storage  since  it  is  necessary 
to  store  coefficients  of  the  least-squares  fits  for  position  angles  from  zero  to  65 
degrees  in  increments  of  five  or  ten  degrees.  Since  these  functions  essentially 
reproduce  the  columns  in  the  range  table,  they  contain  the  deficiencies  inherent  in 
the  range  table;  they  do  not  compensate  for  interaction  among  the  independent 
variables  used  in  the  computation  of  gun  and  fuzing  orders.  Also,  in  the  surface 
mode,  a  non-zero  target  height  is  compensated  for  by  rotating  the  surface  solution 
through  the  position  angle.  This  is  acceptable  for  small  position  angles,  up  to 
approximately  six  degrees,  but  introduces  fairly  large  errors  at  higher  position  angles. 


In  1969,  the  NWL  initiated  a  modest  in-house  research  program  to  further 
investigate  techniques  for  the  implementation  of  gun  ballistics  in  digital  gunfire 
control  computers.  The  program  was  motivated  by  the  anticipated  need  for 
implementing  gun  ballistics  for  long-range  and  guided  projectiles,  a  desire  to 
overcome  the  shortcomings  inherent  in  tae  techniques  used  in  the  MARK  86  GFCS, 
and  the  experience  gained  in  investigating  similar  problems  associated  with  airborne 
digital  fire  control  systems. 


Rather  than  attempt  to  fit  the  columns  of  the  range  table  as  was  done  in 
other  approaches  and  thus  retain  the  deficiencies  inherent  in  the  range  table,  the 
NWL  approach  was  to  first  determine  a  closed-form  approximate  solution  for  the 
basic  equations  of  motion  used  to  mathematically  simulate  the  flight  of  the 
projectile  and  produce  the  range  table.  This  approach  uses  an  “Effective  Drag” 
which  is  determined  for  each  data  point  as  a  function  of  slant  range,  position  angle, 
initial  velocity,  ballistic  density  and  temperature  by  substituting  into  a  20-term 
expression  the  values  of  each  of  these  variables.  This  Effective  Drag  is  then 
substituted  into  the  closed-form  expression  for  elevation  angle.  Time-of-flight  is 
computed  similarly  using  elevation  angle  and  the  variables  mentioned  above. 


After  deriving  the  Effective  Drag  Functions  and  comparing  the  accuracy  and 
storage  with  the  functions  used  in  the  MARK  86  system,  the  Naval  Ordnance 
Systems  Command  (NOSC)  funded  NWL  to  tailor  the  functions  to  5 '754  caliber  and 
3 " / 5 0  caliber  ammunition  for  use  in  the  digital  version  of  the  MARK  68  GFCS,  the 
most  widely  used  naval  gunfire  control  system.  In  addition,  after  discussions  with 
NOSC  and  LEC  personnel,  NWL  and  LEC  were  funded  to  investigate  the  feasibility 
of  replacing  the  ballistic  functions  in  the  MARK  86  system  with  the  Effective  Drag 
Functions.  This  feasibility  study  concluded  that  the  effective  drag  functions  should 
be  used  in  the  MARK  86  system  in  place  of  the  LEC  functions.  Further 
information  concerning  the  feasibility  study  as  well  as  other  significant  conclusions 
of  the  study  can  be  found  in  Reference  7. 
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II.  TRAJECTORY  DATA 


The  trajectory  data  used  as  reference  data  in  the  computation  of  approximating 
functions  were  computed  using  a  two-dimensional  particle  integration  program  (see 
Reference  8  for  an  example  of  this  type  of  program).  The  model  employs  a 
fourth-order  Runge-Kutta  scheme  to  numerically  integrate  the  equations  of  motion 
subject  to  the  following  assumptions'. 

1.  The  earth  is  flat  and  non-rotating  with  gravity  varying  as  a  function  of 
altitude. 

2.  The  projectile  is  a  point  mass. 

3.  Drag  is  the  only  aerodynamic  influence. 

4.  The  atmosphere  is  Navy  Standard  (NAST). 

The  initial  conditions  (initial  velocity,  elevation  angle),  atmospheric  conditions 
(ballistic  density,  ballistic  wind,  ballistic  temperature),  and  terminal  conditions  (flight 
time  for  AA  mode,  target  altitude  for  surface  mode)  were  selected  based  on  the 
conditions  encountered  in  the  tactical  use  of  the  MARK  68  and  MARK  86  systems. 
Table  1  lists  the  type  of  distribution  assumed  for  each  variable  as  well  as  the  range 
of  conditions  used  for  the  functionalization. 


QUO 
3  _  ID 


woow 
Q-<  Qh 

«4-H  <4-1  <4-4 


O  O  O 
O  O  m 
VC'  t  c-- 
-  (N  M 


WWW 

O,  CU  rOn 
m-h  4-h  4-h 


o  o  o 
O  o  in 
co  m  vt 
*— «  (N  (N 


JLj 

3 

C 

.2  ^ 

3 

4-» 

.2 

■M 

>> 

4~> 

ra 

Jh 

4-* 

•  >— 1 

‘o 

O 

C/3 

2 

*00 

C 

5  < 

a 

3 

13 

CD 

> 

3 

M 

> 

PQ 

Q 

Target  Altitude  Uniform  -50/5000  ft 

(Surface  only) 


III.  TRAJ  COMPUTER  MODEL 


To  generate  the  trajectory  data  reflecting  the  distributions  for  each  variable,  a 
computer  model  “TRAJ”  was  developed  by  NWL  personnel.  The  program  is  a 
variation  of  the  basic  particle  trajectory'  integration  program  discussed  in 
Reference  8.  The  only  significant  difference  between  the  models  is  the  method  used 
to  select  initial,  atmospheric,  and  terminal  conditions  for  each  trajectory.  The  basic 
particle  trajectory  model  requires  input  values  lor  each  ol  the  independent  variables 
whereas  the  TRAJ  programs  selects  the  values  for  each  of  these  variables  based  on 
random  sampling  from  the  distributions  assumed  for  each  variable.  To  accomplish 
the  random  sampling,  the  parameters  for  each  of  the  distributions  are  input  to  the 
model  along  with  the  total  number  of  trajectories  to  be  integrated.  After  integrating 
these  trajectories,  the  model  arranges  the  trajectory  data  in  ascending  order  of  slant 
range  (surface)  or  time-of-flight  (AA)  and  prints  out  the  results  of  the  integration 
along  with  the  initial  and  atmospheric  conditions  used  to  generate  the  trajectories. 
In  addition  to  the  usual  trajectory  variables  computed  by  the  model  (slant  range, 
position  angle,  time-of-flight,  and  drift),  the  effective  drag  value  and  the  ratio  of  the 
tangent  of  the  orientation  angle  at  termination  to  the  tangent  of  the  superelevation 
angle  are  printed  out.  The  effective  drag  value  is  computed  using  the  following 

equations:4 


(IJM)2  sin  (2Eg  ) 
XVAC  "  (G)  (Cos  E4*) 


N 


XVAC 

Rh4* 


FK  = 


-  1  -l-  s/3N  -  2 

R4* 


If  a  non-zero  value  of  range  wind  is  generated  by  the  random  sampling 
technique,  the  model  computes  all  of  the  above  data  with  respect  to  the  moving  aii 
mass.  For  further  explanation  of  these  calculations,  see  the  listing  of  the  TRAJ 
model  given  in  the.  appendices.  If  the  target  height  used  to  terminate  the  trajectory 
is  greater  than  the  maximum  ordinate  of  the  trajectory,  the  generated  trajectory  is 
discarded  and  a  new  set  of  conditions  is  generated  using  the  random  sampling 
technique. 


4 A  glossary  of  terms  is  given  in  the  appendices. 
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IV.  BALLISTIC  FUNCTIONS 


The  ballistic  functions  were  derived  from  the  differential  equations  of  motion 
for  a  point  mass.  Since  these  equations  have  no  closed-form  solutions  for 
non-vacuum  trajectories  in  terms  of  elementary  functions  only,  the  solution  is 
obtained  by  approximating  the  true  solution.  The  functions  obtained  by  this 
approximation  are  as  follows  (The  derivation  of  the  superelevation  angle  function  is 
given  in  the  appendices  and  the  derivation  of  the  time-of-flight  function  is  given  in 
Reference  10.): 


V4  =  -sin 

n 


-1 


G  •  N  •  Rh4 


UM2 


T4  =  SQRT 


6  •  A  •  UM  •  sin  Eg*  -  6  •  RV  -  G  •  A2 


2  •  G 


where 


N  =  1  +  (K  •  RS)(2  +  K  •  RS)/3 
Rh4* 


A  = 


UM  •  cos  Eg* 


To  overcome  the  poor  accuracy  of  the  superelevation  equation,  for  long 
times-of-flight,  the  quantity  “K”  is  replaced  by  the  term  “FK”  which  is  expressed 
as  a  least-squares  function  of  initial,  atmospheric,  and  terminal  conditions. 


A  A  mode 


FK  =  K1  +  K2  •  US  +  K3  •  RS  •  US  •  DS  +  K4  •  E4S  +  KS  •  DS  +  K6  ■  US  •  DS 
+  K7  *  RS  •  DS  +  K8  •  US  •  TS  +  K9  •  DS  •  TS  +  K10  •  E4S  •  RS2  + 

Kll  •  RS  •  US  •  E4S  +  K12  •  RS  •  E4S  +  K13  •  US  •  E4S  +  K14  •  RS3  + 

K15  •  E4S3  +K16  •  RS  •  DS  •  E4S  +  K17  •  RS  •  E4S2  +K18  •  US  •  DS  •  E4S  + 
K19  •  RS  •  TS  +  K20  •  E4S  •  US2 


ic  ip 


,1 '  m 


f 


Surface  Mode 


FK  _  K1  +  K2  •  US  +  K3  •  TS  +  K4  ■  E4S  +  K5  •  RS  ■  US  +  K6  •  RS  *  TS 

+  K7  •  US  •  DS  +  K8  *  US  ■  TS  +  K9  •  US  •  E4S  +  K1 0  •  RS2 

+  K11  •  RS  •  E4S2  H-K12  •  US  •  RS2  +K13  •  DS  •  RS2  +K14-  E4S  •  RS2 

+  K15  •  E4S-US2  +  K16  ■  E4S  •  DS2  +  K17  •  RS3  +  K18  *  RS  •  US  •  E4S 

+  K 1 9  •  RS  •  DS  •  E4S  +  K20  •  US  •  DS  •  E4S 

The  accuracy  of  the  time-of- flight  function  is  also  improved  by  the  addition  of  a 

DELT  term  which  is  a  least-squares  function  of  the  initial,  atmospheric,  and  terminal 

conditions. 

AA  Mode 

DELT  =  LI  +  L2  •  E4S  •  RS2  +  L3  *  DS  •  RS2  +  L4  •  US  •  RS2  +  L5  *  RS  •  E4S2 

+  L6  *  RS3  +  L7  •  RS  •  US  +  L8  •  RS  •  TS 


Surface  Mode 

DELT  =  LI  •  US  +  L2  ■  RS  *US  +  L3  •  RS  •  DS  +  L4  •  US  •  TS 

+  L5  •  DS  •  TS  +  L6  •  US  •  RS2  +  L7  ■  DS  •  RS2  +  L8  •  E4S  •  RS2 
+  L9  •  RS3  +L10  •  RS  •  US  •  E4S 

In  addition  to  DELT  and  FK,  two  other  quantities  are  computed  using 
least-squares  expressions.  The  drift  of  the  projectile  due  to  spin  is  computed  using 
the  following  equations: 
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AA  Mode 


DRS  =  Ml  +  M2  •  RS  •  T4S  +  M3  •  RS  •  E4S  +  M4  •  T4S2 
+  M5  •  T4S  •  E4S  +  M6  •  E4S2 


DRIFT  =  DRS  •  10000 


Surface  Mode 


DRS  =  Ml  •  RS  +  M2  •  T4S  +  M3  •  T4S2 
DRIFT  =  DRS  •  10000 


The  fourth  quantity  computed  using  a  least-squares  expression  is  the  ratio  of  the 
tangent  of  the  orientation  angle  of  the  projectile  at  termination  to  the  tangent  of 
the  superelevation  angle.  This  quantity  is  computed  only  in  the  surface  mode  and  is 
used  to  correct  for  Coriolis  effects.  The  following  expression  is  used  to  compute 
this  ratio. 


Ratio  =  N1  +  N2  •  US  •  RS2  +  N3  •  DS  •  RS2  +  N4  •  RS2  •  E4S 
+  N5  ■  RS  •  US  •  DS 


The  cross-terms  in  the  above  expressions  are  used  to  account  for  the  interaction 
of  the  independent  variables.  In  particular,  the  interaction  among  initial  velocity, 
ballistic  density,  and  slant  range  to  target  is  large. 

The  Effective  Drag  Functions  also  have  several  other  advantages  compared  to 
the  MARK  86  ballistic  functions. 


1.  Interpolation  as  a  function  of  target  position  angle  is  not  needed  in  the 
AA  mode  since  the  position  angle  is  one  of  the  independent  variables  used 
in  the  functions. 

2.  Non-standard  values  of  initial  velocity,  ballistic  density,  and  ballistic 
temperature  are  considered  directly  in  the  solution  for  the  superelevation 
angle  and  time-of-flight  rather  than  correcting  the  calculated  values  of  these 
two  quantities  for  standard  conditions. 

3.  In  the  surface  mode,  target  height  is  considered  directly  in  the  solution  for 
the  superelevation  angle  through  the  terms  involving  the  target  position 
angle. 

4.  The  amount  of  storage  needed  for  the  ballistic  functions  and  the  related 
coefficients  is  much  less  than  the  corresponding  storage  for  the  MARK  86 
functions. 

5.  Since  fewer  partitions  of  the  least-squares  functions  are  needed,  the  amount 
of  computer  time  used  to  find  the  correct  partitions  is  less  than  for  the 
MARK  86  system. 

To  compensate  for  horizontal  range  wind,  the  following  procedure  is  used: 

1.  The  horizontal  range  of  the  target  is  modified  using  the  following 

equations 


Rh4*  =  Rh4  -  T4  •  Wrh 


2.  The  initial  velocity  of  the  projectile  is  modified  by  the  apparent  wind 
velocity. 


UM  =  SQRT  (U  sin  Eg)2  +  (U  cos  Eg  -  Wrha)2 


3.  After  computing  a  superelevation  angle  using  the  equation  given  previously, 
the  elevation  angle  is  modified  using  this  equation: 


Eg  =  tan 


UM  ■  sin  Eg* _ 

UM  •  cos  Eg*  +  Wrha 


l 


The  sequence  of  calculations  for  both  the  AA  and  surface  modes  and  a 
glossary  of  terms  are  given  in  the  appendices. 

In  addition  to  the  independent  variables  discussed  above,  many  of  the  variables 
which  effect  the  projectile  trajectory  can  be  taken  into  account  using  closed-form 
expressions  (i.e.,  no  least-squares  coefficients  are  needed).  Lateral  winds  and  the 
affect  of  own  ship  motion  perpendicular  to  the  line-of-fire  on  the  trajectory  can  be 
considered  as  described  in  Reference  11.  Non-standard  projectile  weight  can  be 
considered  by  modifying  the  initial  velocity  of  the  projectile  and  the  ballistic  density 
as  described  in  Reference  1 1 . 


V.  LEAST-SQUARES  PROCEDURES 


The  expressions  previously  given  for  FK,  DELT,  DRIFT,  and  RATIO  were 
determined  by  the  BIVOR  (Backward  Independent  Variable  Ordering  by  Regression 
Sum  of  Squares)  procedure  discussed  in  Reference  12.  This  procedure  ranks 
least-squares  terms  according  to  ascending  order  of  importance.  To  obtain  the 
expression  for  FK,  a  71 -term  least-squares  fit  was  carried  out  using  the  trajectory 
data  generated  by  the  TRA.I  program.  The  least-squares  expression  consisted  of  low 
powers  of  each  of  the  six  independent  variables  and  products  of  these  powers.  By 
deleting  the  least  important  term  at  each  step,  the  twenty  most  important  terms 
were  selected  for  inclusion  in  the  final  functional  form.  The  same  procedure  was 
used  to  select  the  final  functional  forms  for  DELT,  DRIFT,  and  RATIO.  The 
ranking  was  carried  out  using  the  SELECT  computer  model  coded  at  NWL.  The 
major  difference  between  SELECT  and  the  model  discussed  in  Reference  12  is  that 
the  SELECT  program  can  do  least-squares  fits  and  related  analyses  for  several  sets 
of  data  at  one  time.  This  is  accomplished  by  summing  the  values  ol  the  additional 
regression  sum  of  squares  (ARSS)  for  each  term  over  all  the  sets  ol  data  input  to 
the  model.  The  BIVOR  ranking  is  based  on  several  ieast-squares  fits  using  the  same 
functional  form  for  each  data  set.  The  least-squares  coefficients  for  each  data  set  are 
given  in  the  appendices. 

To  attain  the  functional  accuracy  needed  to  adequately  approximate  the 
solution  to  the  particle  equations  of  motion,  it  is  necessary  to  partition  the  fits  as 
a  function  of  slant  range  to  target.  The  AA  functions  required  two  partitions  for 
both  5"/54  and  3"/50  projectiles.  In  the  surface  mode  the  following  number  of 
partitions  are  used: 

5"/54  full  charge  3 

5"/54  reduced  charge  4 

3  "  /  5  0  3 

If  better  accuracy  is  desired,  it  is  possible  to  use  a  larger  number  of  terms  for  each 
of  the  four  variables  fitted.  Terms  different  from  those  currently  used  in  the 
SELECT  or  RANDOM  may  be  added  to  the  programs  with  a  minimum  of  new 
coding. 


FUNCTIONAL  ACCURACY 
AA  MODE 


RMS  in 

RMS  in 

RMS  in 

Gun 

Slant  Range 

Range  Error1 

Miss  Distance2 

Time-of-Flight 

(Yds) 

(Mils) 

(Yds) 

(Sec) 

5"/54 

0-11200 

.8 

.6 

.02 

11200-15000 

.8 

3.0 

.06 

TOTAL 

.8 

1.8 

.04 

3750 

0-8500 

1.8 

2.4 

.04 

8500-11000 

2.1 

7.6 

.10 

TOTAL 

1.9 

4.8 

.07 

OVERALL  TOTAL 

1.5 

3.6 

.05 

’Measured  along  R4 

2Measured  perpendicular  to  R4  at  aim  point 

FUNCTIONAL  ACCURACY 
SURFACE  MODE 


RMS  in 

RMS  in 

Projectile 

Range  Error1 

Time-of-Flight 

(Mils) 

(Sec) 

5754  Full  Charge 

1.9 

.13 

5"/54  Reduced  Charge 

2.5 

.04 

3750 

1.8 

.04 

’Measured  along  R4 


VI.  FUTURE  WORK  IN  BALLISTIC  FUNCTIONALIZATION 

As  discussed  previously,  NWL,  and  LEC  are  currently  tasked  to  modify  the 
MARK  86  tactical  program  tc  allow  the  use  of  file  Effective  Drag  Functions  in 
place  of  the  cubic  polynomials  currently  used.  Since  the  AA  portion  of  the 
MARK  86  program  currently  solves  for  gun  orders  and  fuzing  orders  for  slant  ranges 
exceeding  15,000  yards,  the  effective  drag  functions  will  be  extended  to  include 
slant  ranges  greater  than  the  15,000  yards  limit  now  imposed.  It  also  may  be 
desirable  to  include  more  partitions  for  the  surface  mode  for  the  5'7 54  full  charge 
projectile.  Since  new  5-inch  ammunition  is  currently  being  added  to  the  list  of 
possible  ammunition  for  both  the  MARK  68  and  MARK  86  systems,  it  will  also  be 
necessary  to  fit  the  ballistics  data  for  the  new  rounds. 

In  addition  to  these  modifications,  the  joint  NWL-LEC  study  is  also  aimed  at 
modifying  the  MARK  86  program  to  interface  with  the  8"/55  major  caliber 
lightweight  gun  (MCLWG).  The  trajectory  data  for  the  8-inch  projectile  will  be 
generated  in  the  system  by  use  of  the  Effective  Drag  Functions. 

It  is  envisioned  that  the  Effective  Drag  Functions  can  be  used  in  future  digital 
fire  control  systems  for  projectiles  ranging  in  diameter  from  20  millimeters  to  16 
inches.  There  are  many  advantages  to  having  the  same  functional  form  in  each  of 
the  gunfire  systems:  (1)  Only  one  set  of  computations  is  required  to  generate  the 
least-squares  coefficients  for  each  projectile  type;  (2)  If  ballistics  for  a  round  are 
changed,  new  ballistics  can  be  added  with  a  minimum  of  effort;  (3)  When  analysts 
work  with  the  coding  of  each  of  several  systems,  it  will  not  be  necessary  to  be 
familiar  with  many  different  types  of  ballistics  functions. 

The  anticipated  addition  of  guided  projectiles  to  the  naval  arsenal  of  projectiles 
will  require  ballistic  functions  that  cause  the  projectile  to  pass  through  the 
acquisition  cone  and  thus  acquire  the  target.  The  effective  drag  functions  should 
work  extremely  well  for  this  purpose  since  the  basic  assumptions  used  in  the 
derivation  of  the  effective  drag  functions  are  compatible  with  the  trajectory  data  for 
the  guided  projectile  prior  to  the  actuation  of  guidance. 
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DERIVATION  OF  SUPERELEVATION  ANGLE  EQUATION 


The  particle  equation  for  horizontal  range  is: 


dx  _  d2  x 
dt  ~  dt2 


-KV2  cos  e 


where 


IC  =  >/2  p7CD  is  the  retardation  of  the  particle 
p  is  air  density 

7  is  the  reciprocal  of  the  ballistic  coefficient 
CD  is  the  drag  coefficient 
V  is  the  velocity  of  the  particle 
0  is  the  orientation  angle  of  the  particle. 

Substituting  ds/dt  for  V  gives 


dx 

dt 


-KV  cos  0 


ds 

dt 


or 


dx  =  -  KV  cos  6  ds. 


Substituting  x  for  V  cos  0  in  (3)  gives 


dx  =  -Kxds. 


Separating  variables  gives: 


dx 

x 


=  -  Kds . 


A-l 


Integrating  (5)  and  substituting  limits  gives 


log(V  cos  <p)  -  log  x  =  ICs 


where 


(p  is  the  initial  value  of  6. 


log(V  cos  (p/x)  =  ICs , 


V  cos  (p/x  =  eKs  , 


and 


=  V  cos  (pe~KS  =  xne~ 


Ks 


Also 


V2 

—  -  g  cos  6 


where 


r  is  the  radius  of  curvature,  and 
g  is  the  acceleration  due  to  gravity. 


Since  V  =  x  sec  6  and  r  = 


-ds 

dfl"’ 


equation  (10)  becomes 


x2  sec2  d 
-ds/d  0 


g  cos  0 


or 


x2  sec2  0  =  (~g  cos  d) 


ds 

dd' 


Therefore 


ngpim 


RPMH 


sec2  6  d 6  - 


■  g  cos  0  ds 


Substituting  dx  for  ds  cos  0  gives: 


sec2  Odd  - 


Substituting  the  expression  from  (9)  into  Equation  (14)  gives: 


sec20d0  - 


V2  cos2  0e  2K  s 


sec20d0 


-  ge2  K  s  dx 
V2  cos2  0 


Using  the  approximation  s  -  x  in  (16)  gives. 


sec20d0  = 


-ge2Kxdx 
V2  cos2  0 


Integrating  (17)  between  appropriate  limits  gives: 


tan  0  -  tan 0  = 


_ 1 _  e2Kx  -  1 

_2KV2  cos2  0 J  l 


(2Kx)2  (2Kx)3 

e2Kx  =  1  +  2Kx  +  — -  +  ~ — 


e2Kx  _  i +2Kx  +  2K2x2  +  -K3x3 


,2Kx  _  ,  _  2Kx  +  2K2 x2  +  -K3x3 


iii a , ii  l  up  ti  1‘iijiatiii  ;r  i  i 1  ft  s  im  m  ivurm ;  *k  io  md 


vv  !'.»v:  1 .  'H  If ’r^nlrflf!-! W Hr :T : !l« 


“f. 


e2Kx  -  1  =  2Kx  1  +  Kx  +  - K2x2 


Substituting  this  expression  into  (18)  gives: 

dy  .  „  f  g 


-  tan  6  =  tan  0 


g 

~  ~ivTn - 7T  2Kx(l 

_2KV~  cos2  0j 


+  Kx  +  -  K2  x2 ) 


gx  9 

dy  =  *an0-  — - —  (, +Kx+  -K2x2)  dx. 

L  V-  cos2  0  3  2 


Integrating  (21)  gives: 


y  =  x  tan  0  ■ 


gx2 

r  2Kx 

K2x2' 

2V2  cos2  0 

1  _l_  _ 

3 

+  - 

3 

Let  N  = 


2Kx  K2x2 

+  -  + - 

3  3  ' 


-et  y  =  0: 


y  =  x  tan  0  ■ 


2V2  cos2  0 


x  tan  0  =  — ^ - 

2V2  cos2  0 

(2V2  cos2  0)(x  tan  0)  =  Ngx2 
2V2  cos2  0  tan  0  =  Ngx 
2V2  cos  0  sin  0  =  Ngx 


‘►1  *«WM 


sin  2 <p 


Ngx 
V2  ' 


(27) 


Therefore 


APPENDIX  B 

SEQUENCE  OF  CALCULATIONS 
AA  AND  SURFACE  MODES 


HwitSIwipiS-w 


pSl*! 


Hi 

fllillflfillFffllffllfBltlll'!'' 


SEQUENCE  OF  CALCULATIONS 


AA  MODE 


Rh4*  =  Rh4  -  T4  •  Wrh 


(1) 


R4*  =  SQRT[(Rh4*)2  +  (Rv)2  ] 


(2) 


E4*  =  tan  1 


(3) 


UM  =  SQRT[(U  sin  Eg)2  +  (U  cos  Eg-  Wrha)2  ]  (4) 


(Use  best  estimates  of  T4  and  Eg  available  at  this  time.  If  no  estimates  are 
available,  assume  Wrha  =  0.) 


FK  =  K1  +  K2  •  US  +  K3  •  RS  •  US  •  DS  +  K4  ■  E4S  +  K5  •  DS  +  K6  •  US  •  DS 
+K7  •  RS  •  DS  +  K8  •  US  •  TS  +  K9  •  DS  •  TS  +  K10  •  E4S  •  RS2  + 

IC1 1  •  RS  •  US  •  E4S  +  K12  •  RS  •  E4S  +  K13  •  US  •  E4S  +  K14-  RS3  + 

Id  5  •  E4S3  +  K16  •  RS  •  DS  •  E4S  +  K17  •  RS  •  E4S2  +  K18  •  US  •  DS  •  E4S  + 
K19  •  RS  •  TS  +  K20  •  E4S  •  US2  (5) 


N  =  1  +  (FK  •  RS)(2  +  FK  •  RS)/3 


(6) 


V4*  =  1/2  sin  1 


N  •  Rh4*~ 
UM2 


(7) 


Eg*  =  E4*  +  V4* 


(8) 


DELT  =  LI  +  L2  ■  E4S  •  RS2  +  L3  •  DS  •  RS2  +  L4  •  US  •  RS2  +  L5  ■  RS  •  E4S2 
+  L 6  •  RS3  +  L7  •  RS  •  US  +  L8  •  RS  •  TS  (10) 


T4  =  SQRT 


-  Rh4* 

UM  •  cos  Eg* 

A  •  UM  •  sin  Eg*  -  6  •  Rv  -  G  •  A2 
2  ■  G 


+  DELT 


(ID 


(12) 


i 

! 


SEQUENCE  OF  CALCULATIONS 
SURFACE  MODE 


=  N1  +  N2  •  US  •  RS2  +  N3  •  DS  •  RS2  +N4  •  E4S  •  RS2  +  N5  •  RS  •  US  •  DS  (1) 


DELX  = 


TANW  =  Ratio  •  tan  Eg 


18274  -  US2  •  cos  (2  •  Eg) 

3  +  4  ■  FK  •  RS  +  3  •  FK2  •  RS2 


+  D1 


AG  -  [0.2431  •  T4]  [.1038  •  DELX  +  1.8  •  RS/TANW] 


BG  -  [0.3646  •  T4  •  RS] 


3  •  tan  Eg  +  TANW 
[_  tan  Eg  +  TANW  J 


DG  =  [0.1215  •  T4  •  RS  •  tan  Eg] 


19  •  tan  Eg  +  TANW 


7  •  tan  Eg +3  •  TANWJ 


DELXO  =  AG  •  cos  L  -  sin  Az 


DOMEGA  =  BG  •  sin  L  -  DG  •  cos  L  cos  Az 


Rh4*  =  Rh4  -  DELXO  -  T4  •  Wrh 


R4*  =  SQRT[(Rh4*)2  +(Rv)2] 


E4*  =  tan 


,-i 


Rv 


Rh4* 


Um  =  SQRT[(U  sin  Eg)2  +  (U  cos  Eg  -  Wrha)2 


(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(ID 

(12) 


(Use  best  estimates  of  T4,  Eg,  and  FK  available  in  (2),  (3),  (5),  (6),  (9),  and 
(12).  If  none  are  available,  assume  Wrh,  DELXO,  and  DOMEGA  are  zero.) 


=  K1  +  K2  •  US  +  K3  *  TS  +  K4  •  E4S  +  K5  •  RS  •  US  +  K6  •  RS  •  TS 


+  K7  •  US  •  DS  +  K8  •  US  •  TS  +  K9  •  US  •  E4S  +  K10  ■  RS2 


+  K1 1  •  RS  •  E4S2  +  K12  •  US  *  RS2  +  K13  •  DS  •  RS2  +  K14  •  E4S  •  RS2 


+  K15  •  E4S  •  US2  +  K16  •  E4S  •  DS2  +  K17  •  RS3  +  K18  •  RS  •  US  •  E4S 


+  K 1 9  •  RS  •  DS  •  E4S  +  K20  •  US  •  DS  •  E4S 


N  =  1  +  (FK  •  RS)(2  +  FK  •  RS)/3 


sin  2  = 


G  •  N  •  Rh4* 


(If  sin2  is  greater  than  one,  limit  it  to  one.) 


V4*  =  1  /2  sin" 1  (sin2) 


(If  a  high  angle  solution  is  desired,  set  V4*  -  1/2[180°-  sin  1  (sin2)]). 


Eg*  =  E4*  +  V4* 


Eg  =  tan" 


Um  •  sin  Eg* 

Uni  •  cos  Eg*  +  Wrha 


DELT  =  LI  •  US  +  L2  •  RS  •  US  +  L3  •  RS  •  DS  +  L4  •  US  •  TS 


+  L5  •  DS  •  TS  -I-  L6  •  US  ■  RS2  +  L7  *  DS  •  RS2  +  L8  •  E4S  •  RS2  (19) 


+  L9  •  RS3  +  L10  •  RS  •  US  •  E4S 


SWS 


APPENDIX  C 
GLOSSARY  OF  TERMS 


GLOSSARY  OF  TERMS 


A 

AG 

AZ 

BG 

DELT 

DELX 

DELXO 

DG 

DOMEGA 

DRIFT 

DRS 

DS 

D1 

Eg 

Eg* 

E4* 

E4S 

FK 

G 

K1  -  K20 
L 

LI  -  L10 


Intermediate  quantity  used  in  computation  of  T4 
Intermediate  Coriolis  coefficient 

Azimuth  of  the  line  of  fire  (measured  clockwise  from  the  north) 

Intermediate  Coriolis  coefficient 

Intermediate  quantity  used  in  computation  of  T4 

Change  in  range  corresponding  to  an  increase  of  ten  minutes  in 
Eg  (Yards) 

Change  in  range  due  to  Coriolis  effect  (Feet) 

Intermediate  Coriolis  coefficient 
Cross  range  due  to  Coriolis  effect  (Feet) 

Drift  (Feet) 

Drift  X  10~ 4  (Feet) 

Ballistic  density  X  10" 2  (Percent  of  Standard) 

Coefficient  in  DELX  expression 

Elevation  angle  with  respect  to  ground  (Degrees) 

Elevation  angle  with  respect  to  air  (Degrees) 

Aiming  position  angle  (Degrees) 

Aiming  position  angle  (Radians) 

Intermediate  quantity  used  in  computation  of  Eg 
Gravity  constant  (32.174  Ft/Sec2) 

Coefficients  for  FK .  expression 
Latitude 

Coefficients  for  DELT  expression 


C-l 


0ptriS*W  ?>«****“ 


Ml  —  M5 


APPENDIX  D 


COEFFICIENT  TABLES 


TABLE  OF  COEFFICIENTS 


A  A  MODE 


5754 


3750 


0-11200  yards  11200-15000  yards  0-8500  yards 
Slant  Range  Slant  Range  Slant  Range 


8500-11000  yards 
Slant  Range 


K1 

-.9810293E-01 

K2 

.2544428E-00 

K3 

-.2724356E-01 

K4 

.4287342E-00 

K5 

.4066069E-00 

K6 

-.7607970E-00 

K7 

.1738324E-01 

co 

^4 

-  .1 100654E-01 

K9 

.1377282E-01 

K10 

.2904165E-02 

Id  1 

-.2465604E-00 

K12 

.1 197910E-01 

K13 

-. 277737 1E+01 

K14 

.3755598E-04 

IC15 

-.9235122E-02 

K16 

.2838094E-01 

K17 

.7251661E-02 

K18 

-.3073384E-00 

K19 

.1269905E-02 

K20 

.6224520E+01 

LI 

-.9330092E-02 

L2 

-.5993651E-01 

L3 

.4183774E-01 

IA 

.3654574E-02 

L5 

.4870824E-01 

L6 

;6214550E-03 

L7 

-.8304742E-01 

L8 

.3393555E-02 

-.1874708E-00 
.58777 15E+00 
.1  1 53756E-00 
.5430970E+00 
.5890901 E+00 
-.1  574866E+01 
-.7027828E-02 
.4833599E-01 
.79298 18E-02 
.1208286E-01 
-.845 1241  E+00 
.6633662E-01 
-.4274047E+01 
-.2408244E-03 
-.225421 2E-  01 
.5660594E-01 
.1474725E-01 

-  .701 8579E+00 
-.8923626E-03 

.1371892E+02 

.3839703E-00 
-.101 1252E-00 
-.85093 10E-02 
.8007668E+00 
.1231386E-00 

—  .1955255E— 01 
-.1392877E+01 
-.1643145E-01 


.4234776E-00 

-.2170355E+01 

.8540107E+00 

-.2514494E-00 

.7784072E+00 

-.1034142E+01 

-  .9297848E-01 
.41 42487 E- 00 

—  .47 1 597 1 E-01 
.37 1 4500E-  0 1 

-.1279363E+01 
.1255800E-00 
.1783080E+01 
-.4135671E-02 
-.3473544E-02 
.9366087E-01 
.7283655E-02 
-.3798729E-00 
— .  1 24581 6E-0 1 
-.5342851E+00 

.242697  IE- 00 
-.1 199946E-00 
-  .521 5645E-01 
.3709750E+01 
.659091 5E- 01 
-.2250537E-00 
-.3104506E+01 
-.5360088E-01 


.3488444  E- 00 

-  .1 5853 19E+01 
.703  249  5  E+00 

-.7887737E+00 

.8533952E+00 

-  .6907776E+00 
-.1725396E-00 

.7856080E+00 
-.801471  IE- 01 
.1300060E-00 
-.4156714E+01 
-.1352547E-00 
.7640870E+01 

-  .1 502046E-03 

-  .5204026E-01 
.6  8 19491  E+00 
.5074930E-01 

-.525541 1E+01 
-.3863329E-01 
.9673920E+01 

-.1428936E+01 
-.2019259E-00 
-.2341028E-00 
-.3098614E-00 
.164261  IE-00 
-  .1229 106E-01 
.6509672E+01 
-. 74403 88E- 01 


Ml 

M2 


.348241 8E-03  .3772293E-02 

—  .1 47 1 266E— 01  -.2745239E-01 


.4170855E-04  .6598455E-02 

-.2647575E-  01  -  .7894141E-01 


TABLE  OF  COEFFICIENTS  (Continued) 

5754  3750 

0-11200  yards  11200-15000  yards  0-8500  yards  8500-11000  yards 
Slant  Range  Slant  Range  Slant  Range  Slant  Range 

.2738276E-02  .5895750E-02  .6164067E-02  .1671360E-01 

.4269650E-00  .5473077E+00  .6476893E+00  .1 122133E+01 

-.6064290E-01  -.9502992E-01  -.9995276E-01  -.1934456E-00 

-.3310132E-03  -.3940700E-02  -.1 100139E-02  -.7420265E-02 


TABLE  OF  COEFFICIENTS 


SURFACE  MODE 


5"/54  FULL  CHARGE 


0-18000  yards 

18000-24000  yards 

14500-24000  yards* 

Slant  Range 

Slant  Range 

Slant  Range 

K1 

.3955 169E-00 

.5669702E+00 

-.1704694E-00 

K2 

-.191 2286E+0 1 

-.3534244E+01 

.4507427E+01 

K3 

-.2172905E-01 

-.3670860E-01 

-.2749037E-01 

K4 

-.2803963E+01 

-.8572506E+01 

-.8533351 E+00 

K5 

-. 3226089 E-01 

.2620820E-  00 

-.121 734 1  E+0 1 

K6 

-.2404475E-02 

—  .5 1 3 1 253E-02 

-.1922547E-02 

K7 

.1 08982 1E+01 

.1 543458E+01 

.4899658E-00 

K8 

.1 749775E-00 

.2731 868E-00 

.1  530299E-00 

Iv9 

.2349699E+02 

.6977835E+02 

.1272636E+02 

IC10 

.237191  IE-02 

-.4980339E-02 

.1 525784E-01 

IC]  1 

.2128974E-01 

.4630666E-01 

-.1 664827E-00 

K12 

.2455250E-01 

.1 387262E-01 

.7428142E-01 

K13 

.2164224E-03 

-  .4995185E-02 

.1 606694E-02 

K14 

.5628305E-02 

.204201 5E-01 

-.4869881E-02 

K15 

-.2933482E+02 

-.1  508217E+03 

-.7734495E+02 

K16 

.1004287E+01 

—  .135461 9E+0 1 

.1  123986E-00 

K17 

-.1071262E-02 

.1 1 21 990E-04 

-.1 394350E-02 

K1 8 

-.34 1636  IE- 00 

-.1  131332E+01 

.2505003E+01 

Id  9 

.44303 17E-01 

.826351 0E—  0 1 

-.8141508E+00 

K20 

-.8659542E+01 

.969141 5E+01 

.1 692769E+02 

LI 

.7572007 E+00 

.9917354E+01 

.3534241 E+02 

L2 

-.1071067E+02 

-.9947143E+01 

-.3416829E+01 

L3 

.22503 18E+01 

.1  541953E+01 

-.381 25 1 5E+01 

L4 

.7172221 E+0 1 

-.1404480E+02 

-.6 125264E+01 

L5 

-.1951473E+01 

.21251 57E+01 

-.6025834E-01 

L6 

.2888127E+01 

.2  886421  E+0 1 

-.5432146E+00 

L7 

-.4459809E-00 

-.4308659E-00 

.6589226E+00 

L8 

-.9974401 E-01 

-.9333276E+00 

.4681561E-00 

L9 

-.3265210E-01 

-.2801638E-0  1 

-.58101 61E-02 

L10 

.51 16443E-01 

.1 800316E+02 

-.1704057E+02 

High-anglc  partition 


D-3 


TABLE  OF  COEFFICIENTS  (Continued) 


0-18000  yards 

18000-24000  yards 

14500-24000  yards 

Slant  Range 

Slant  Range 

Slant  Range 

11 

2657291 E-02 

—  .1615951 E—  0 1 

-.21 56430E-01 

12 

.41 56908E-01 

.2532693E-00 

.3528588E-00 

13 

.1 81 5574E-00 

.9 500279 E- 01 

.2902570E-01 

11 

.9191 189E+00 

.  1340*61 E+0 1 

.1 1 33630E+01 

12 

-.2446265E-01 

.9275399E-01 

.5260185E-01 

13 

-.3669022E-02 

-  .44726 10E-  01 

-.5129406E-01 

14 

.1 9871 27E-01 

.5253587E-01 

.1678999E+00 

15 

.1005881 E+0 1 

.8647075E+00 

.1342037E+01 

High-angle  partition 


TABLE  OF  COEFFICIENTS 
SURFACE  MODE 


5"/54  REDUCED  CHARGE 


0-7000  yards 

7000-9000  yards 

9000-13500  yards 

Slant  Range 

Slant  Range 

Slant  Range 

ICl 

.5298128E+00 

.12021 59E-  00 

.2826655E-00 

K2 

-.3934450E+01 

-.3309993E-00 

-.5596140E+01 

K3 

-.1  561971E-00 

-.1674386E-00 

-.9466058E-01 

K4 

-.1797008E+02 

-.8886847E+01 

-.2815730E+01 

K5 

.3823152E-00 

-.1078237E-01 

.4734604E+01 

K6 

-.2544909E-01 

-. 63  209  5  3  E- 02 

.2677492E-03 

K7 

.1 948643 E+01 

1780902E+01 

.6738246E+00 

K8 

.1264884E+01 

.1072910E+01 

.4546266E-00 

K9 

.2505323E+03 

.1281992E+03 

.4029905E+02 

K10 

-.8957101E-01 

-.1 1 69 131 E—  0 1 

-.1 126286E-00 

Kll 

-.1008650E-00 

.3548671E-01 

.2144906E-00 

K12 

.582081 6E+00 

-.3706394E-02 

-.7773864E+00 

K13 

-.1608779E-01 

-.12381 16E-01 

.3467924E-02 

ICl  4 

-. 308680  IE-02 

-.249499  IE- 00 

.3212877E-00 

ICl  5 

-.8940509E+03 

-.5395590E+03 

-.68521 14E+02 

K16 

.3083025E-00 

-.7471392E+00 

-.468205  IE- 00 

K17 

.2159488E-04 

.5139764E-02 

.2243239E-01 

K18 

.4008211 E+01 

.79  8447  5  E+01 

-.1207549E+02 

K19 

-.5602425E+00 

.572871  IE-01 

-.9943378E-01 

K20 

.2061799E+01 

.7491 158E+01 

.1 1 12658E+02 

LI 

.781 1849E+00 

.1045614E+02 

-.1795124E+02 

L2 

-.2800263E+01 

.5766946E+01 

.8842858E+01 

L3 

.1389815E-00 

-.2489663E+01 

.1 581421E-00 

L4 

44971 83 E- 00 

-.7774090E+01 

-.6966992E+01 

L5 

-.4123882E-02 

.1035831E+01 

.1233712E+01 

L6 

.4628736E+01 

-.1 142702E+01 

.2951062E-00 

L7 

-.3203420E-00 

.8058107E+00 

-.2294523E-00 

L8 

-.34 1252  IE- 00 

.1092594E-00 

-.7223061 E+00 

L9 

-.1430581E-00 

-.1 131578E-00 

—  .3 1 85 181 E-  01 

L10 

.1875993E+01 

-.4491435E+01 

.1073504E+02 

Ml 

-.3533008E-02 

-.6760702E-02 

-.1267627E-01 

M2 

.4561825E-01 

.7988161E-01 

.1492779E-00 

M3 

,190720  IE- 00 

.1994128E-00 

.1 620066E-00 

N1 

.9594301 E+00 

.9225252E+00 

.1050865E+01 

N2 

.5544641 E+00 

.1 160196E+00 

.49231 14E- 01 

N3 

-.1 597535E+00 

-.7929263E-01 

-.2162164E-01 

N4 

.1302047E+00 

.1441210E+00 

.1523246E+00 

N5 

.2194533E+01 

.204793  5E+01 

.8052037E+00 

*High-angle  partition 


D-5 


7000-13500  yards* 
Slant  Range 

-.7421792E+00 
.1 540729E+02 
-.6586182E-01 
-.3939762E-00 
-.7041342E+01 
.751  1573E-02 
.6031 804E+00 
.1712773E-00 
.2441291E+02 
.128007  IE- 00 
-.5780454E+00 
.924239  5  E+00 
.5279097E-02 
.4976959E-01 
-.1818477E+03 
.364694  IE-00 
-.2476135E-01 
.44266 19E+01 
-.1679682E+01 
.27525  67E+02 

-. 294760 1E+02 
.3014125E+01 
.1 148676E+01 
-.1046813E+02 
.1 567005E+01 
-.3058953E+01 
-.2975420E-00 
.82368 17E+00 
.124961 5E- 00 
-.1541037E+02 

-.1 881960E-01 
.2039353E-00 
.1337699E-00 

.1454356E+01 
.1853321E+00 
-.1884167E+00 
.5649786E+00 
.3291 107E+01 


TABLE  OF  COEFFICIENTS 


SURFACE  MODE 


3"/50 


0-7500  yards 

7500-14200  yards 

Slant  Range 

Slant  Range 

K1 

.9972894E+00 

.1269381  E+01 

K2 

-.4329356E+01 

-  .7409070E+01 

K3 

-6021274E-01 

-.1221252E-00 

K4 

-.1048994E+02 

-.9322561E+01 

K5 

-.1499602E-01 

.2043580E+01 

K6 

-.9419808E-02 

—  .1481 877E— 0 1 

K7 

.2133987E+01 

.2687206E+01 

K8 

.4405334E-00 

.6623383E+00 

K9 

.8209066E+02 

.70059 17E+02 

K10 

-.5925505E-02 

-.1004642E-00 

Kll 

-.5723983E-01 

.5757654E-01 

K12 

.1 540559E-00 

-.2100219E-00 

K13 

.5740481 E- 01 

.9732245E-03 

K14 

-.1439001E-01 

.1  546526E-00 

K15 

-.1 566121E+03 

-.1 127932E+03 

K16 

.44721 81 E+01 

-.3886304E-00 

K17 

-2340028E-01 

.1453335E-01 

K18 

.8741858E+01 

-.4501 337E+01 

K19 

-.2226840E+01 

.3944278E-00 

K20 

-.1827209E+02 

.3684374E-00 

LI 

.5494928E+00 

.3971734E+01 

L2 

-.6203279E+01 

-.7236301 E+00 

L3 

.1107009 E+01 

-.6270852E+00 

L4 

.1432898E+01 

-.2510048E+01 

L5 

-.4216346E-00 

.426392 IE- 00 

L6 

.4167495E+01 

:1042376E+01 

L7 

-.4582454E-00 

.1 176264E-00 

L8 

-.4837642E-00 

-.7669436E+00 

L9 

-.1492970E-00 

-.5601244E-01 

L10 

.3147546E+01 

.6247304E+01 

^High-angle  partition 

D-6 


It 

■ 

•j 

9000-14200  yards* 

Slant  Range 

-.48751 50E-00 
.1 34599 1E+02 
-.5220929E-01 
-.4927079E+01 
-.7763816E+01 
.21 58699E-02 
.3029355E+01 
.1581798E-00 
.5456838E+02 
.2144224E-00 
-.5734877E+00 
.1047363E+01 
-.2429330E-01 
-.2474639E-01 
-.1460537E+03 
-.1067673E+01 
-.3850347E-01 
.3077926E+01 
- .  1 607  5 1 3E+0 1 
.2229303E+  02 

.  146361 3  E+02 
-.1614197E+02 
-.1 149805E+01 
-. 3478743  E+01 
.67071 73E+00 
.2648430E+01 
.8012168E+00 
.3390402E-00 
-.1008407E-00 
-.9377208E+01 


TABLE  OF  COEFFICIENTS  (Continued) 


0-7500  yards 

Slant  Range 

7500-14200  yards 

Slant  Range 

9000-14200  yards' 
Slant  Range 

Ml 

M2 

M3 

-.2179687E-03 

-.2067184E-02 

.2951480E-00 

-.920181  IE-02 
.1073744E-00 
.3104029E-00 

-.3402695E-01 
.403 1 741 E- 00 
.99209  5  5  E- 01 

N1 

N2 

N3 

N4 

N5 

.8736376E+00 

-.8522887E-01 

-.3916639E-02 

.1339467E+00 

.2333159E+01 

.1 255767E+01 
.9501215E-01 

76645 64E- 01 
.1909898Er00 
.1736291E+01 

.1353942E+01 
.3830326E+00 
-.2977684E+00 
.93985 17E+00 
.3498572E+01 

*High-angle  partition 
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